Ring-shaped luminescence pattern in biased quantum wells studied as a steady state 

reaction front 
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Under certain conditions, focused laser excitation in semiconductor quantum well structures can 
lead to charge separation and a circular reaction front, which is visible as a ring-shaped photolu- 
minescence pattern. The diffusion-reaction equations governing the system are studied here with 
the aim of a detailed understanding of the steady state. The qualitative asymmetry in the sources 
for the two carriers is found to lead to unusual effects which dramatically affect the steady-state 
configuration. Analytic expressions are derived for carrier distributions and interface position for a 
number of cases. These are compared with steady-state information obtained from simulations of 
the diffusion-reaction equations. 



I. INTRODUCTION 



In mid-2002, two semiconductor-optics experimen- 
tal groups reported dramatic ring-shaped photolumines- 
cence patterns when a focused laser was used to excite 
electron-hole pairs near a coupled quantum well system 
biased with an electric field [1, 2]. Despite initial specu- 
lation invoking Bose-Einstein condensation of excitons, 
it was later found that the luminescence ring is well- 
explained by classical reaction-diffusion dynamics of the 
electrons and holes [3, 4]. The idea is that there is net 
hole injection into the quantum well near the laser ir- 
radiation spot, together with an electron source due to 
leakage current that is roughly uniform across the two- 
dimensional (2D) quantum well plane. This combination 
can lead to a charge-separated steady state configuration, 
with a circular hole-rich island sustained by the localized 
hole source in an electron-rich sea. The interface, where 
outward-diffusing holes recombine with inward-diffusing 
electrons, is the luminescence ring. 

The position of the interface, i.e., the radius of the 
luminescence ring, is not well-understood theoretically, 
despite some theoretical [4, 5] and experimental [5, 6] ef- 
forts. While a full understanding may or may not require 
extra ingredients in addition to the diffusion-reaction 
model [5], a thorough study of the behavior of the in- 
terface position within the diffusion model is certainly 
a necessary first step. The present Article fills this gap 
by presenting a detailed analysis of the steady state, ad- 
dressing aspects such as the position and width of inter- 
face, density distributions, etc. There are a number of 
length scales in the problem which we identify cleanly. 
The phenomenon is put into the context of previous the- 
oretical studies of steady-state reaction fronts and varia- 
tions thereof [7-15]. By considering various possible rela- 
tive values of the tunneling decay rates of the two carrier 
species, we clarify the roles of the tunneling strengths 
in determining the steady-state configuration. A curious 
feature of the steady state is that the reaction zone has, 
in addition to the sharp interface, an extended feature 
on one side where the luminescence does not vanish but 



instead is a nonzero constant. This aspect turns out to 
have a drastic influence on the interface position and the 
overall steady-state structure, which we explain in detail. 

In comparison with previous theoretically studied 
diffusion- annihilation systems [7-15], the present prob- 
lem has several unusual features which justify an ex- 
tended study. These include the single-particle (tunnel- 
ing) decay of one or both species, and the fact that one 
of the reacting species has a spatially extended source 
spanning both sides of the interface. In addition, while 
diffusion-controlled reaction interfaces and patterns have 
been studied in a wide variety of chemical, biological and 
fluid flow contexts [16-18], they are rather novel in elec- 
tronic systems. Indeed, this may well be the only cur- 
rently known example of a diffusion-limited nonequilib- 
rium reaction front or pattern in electronic systems. Fur- 
thermore, there is the intriguing possibility of studying 
quantum phenomena in the ring region [19, 20], where 
the carriers have had time to cool down to quantum de- 
generacy. 

Sec. II introduces the diffusion-reaction equations, sim- 
plifying source details, and also presents the important 
length scales. Sec. Ill contains the analysis of the steady 
state and the main results of this Article. In Sec. IV, 
we point out some limitations of the current model and 
put the current project into context by briefly review- 
ing the relevant theoretical literature. In Sec. V, our 
calculations on the ring radius are put into perspective 
by discussing experimental issues and other calculations. 
The method used for numerical evolution of the diffusion- 
reaction equations is outlined in App. A. 



II. DIFFUSION-REACTION MODEL 

For experimental details beyond what is sketched here, 
the reader is referred to Refs. [1, 2, 6, 22]. The phe- 
nomenon occurs in a two-dimensional quantum well sys- 
tem, either a single well or two closely separated parallel 
wells. Electron-hole pairs are created in the vicinity of 
the well(s), mainly in the substrate, by a focused laser 
excitation. 
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FIG. 1: Top panel: steady state carrier density distributions, 
rih dominates at small r and rie dominates at large r. Bottom 
panel: luminescence profile. Dashed and solid lines corre- 
spond to Pe 7^ and Po = respectively. The hole excess 
Ph — Pe is the same in the two cases, so that the only dif- 
ference is in the structure around r = 0. Lower-right inset: 
band-structure schematic (single well). Arrows indicate tun- 
neling processes described in text. 

A voltage is applied across the well(s) using conducting 
electron-rich (n+) regions on both sides of the well(s) as 
leads. The original motivation was to enhance the life- 
time of excitons or electron-hole gases by spatially sepa- 
rating electrons and holes in the direction transverse to 
the well(s). A band-structure cartoon of the experimen- 
tal setup is shown in the inset of Fig. 1. Due to the 
electric field bias, there is an influx current of electrons 
into the well as well as a tunneling-out process. In addi- 
tion, the holes can tunnel out in the other direction; this 
corresponds to an electron from the left lead or substrate 
filling up one of the hole states in the well. The three 
processes are shown by arrows in the Fig. 1 inset. Note 
that there is no source of holes due to the biasing. Holes 
are only created by photo-excitation. 

Incorporating the above effects, one can write down 
two-dimensional diffusion-annihilation equations for the 
densities of holes (rih) and electrons (rio) within the quan- 
tum well(s): 

dflii o _„2 /j2 1 

— - = I^hV^rih + Phe - juhno Uh 

dt Th 

(la) 

111 1 

— ^ = DcV^nc + P^e-" - 7nhnc -f G 

at To 

(lb) 

The i?h,o are diffusion constants. The nh,e/''h,o terms 
model decay due to carriers tunneling out of the well(s), 



the r's being tunneling lifetimes. G is the spatially uni- 
form source term for electrons, which is absent for holes. 
The 7 terms represent electron-hole recombination. The 
Ph,G terms are the laser excitation terms; the carriers are 
focused onto a spot roughly of radius ^l- 

By numerically evolving Eqs. (1) in time, one can de- 
termine the steady-state carrier distributions and lumi- 
nescence after Ph.c are turned on. A typical steady- 
state distribution, obtained with > Pc, is shown in 
Fig. 1. The simulation (App. A) is one-dimensional, so 
that radial plots are sufficient. The steady state displays 
a species-separated configuration, together with a peak 
in the luminescence marking the interface, as described 
previously. For Po 7^ 0, the luminescence profile also 
shows an inner peak near r = 0, corresponding roughly 
to a central luminescence spot observed experimentally. 

Typical values of the parameters are taken to be of the 
following orders, D's: several cm^/s, r's: 10~^ s, 7: 10"'^ 
s-icm2, Zl: lO'^ cms, G: 10^^ s-^cm-^, and PZ^'s: lO^^ 
s~^. Units will be omitted in the rest of this Article. 

Eqs. (1), with several variations, was proposed in 
Refs. [3, 4] as the luminescence ring mechanism, and 
studied further in Refs. [5, 6, 19-21]. For a restricted 
case, Ref. [4] also contains a minimal analytic treatment 
of the steady state. 

For the charge separation phenomenon, we need more 
holes diffusing out of the excitation region than electrons. 
In previous studies [3-6], the philosophy has been to in- 
voke differences of unknown origin in the efficiency of 
accumulation in the well(s), i.e., to use Ph > Pe with- 
out detailed explanation. The current understanding of 
carrier asymmetry is thus unsatisfactory. In fact, it is 
possible to have an excess of holes and a resulting lumi- 
nescence ring with Ph = Pq. However, the present author 
will postpone to a future publication an analysis of the 
source asymmetry and of the inner spot structure. 

Since we neglect the inner structure in this study, 
it is convenient to drop the electron source altogether 
(Pc — 0), and assume a point source for the holes, i.e., 
Phe^'' /'l is replaced by P^8(v). For comparison with 
the numeric simulations, where a finite ?l has been used, 
the correspondence is Px = Tr^L-Ph- This is obtained by 
equating outward flux for the point and gaussian sources. 
One result of omitting the electron source is the ab- 
sence of an inner luminescence spot (Fig. 1, solid line 
in lower panel). Moreover, the expressions for hole den- 
sity in Sec. HI will diverge at the illumination spot. This 
(minor) un-physical result is a result of the un-physical 
"point" source. 

We now identify the length scales present in the prob- 
lem. The two most important ones are the depletion 
lengths for electrons and holes, = VDcTc and Zh = 
V-DhTh. The depletion lengths provide the length scales 
for the variation of steady state densities, analogous 
to the diffusion lengths \/Di in the literature on time- 
dependent front formation between two initially sepa- 
rated reactants [8-11], where \/Di gives the spatial vari- 
ation length scale after time t. 
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The ratio of the source strengths, and G, provide s 
a third length scale, which we define as Zgi-c = \/ -Fx/ttG . 
The radius of the ring-shaped interface increases mono- 
tonically with the length 4rc- The interface radius it- 
self, and the interface width are not input parameters 
in the problem but emerge from the analysis as impor- 
tant length scales. We are interested in cases where the 
interface is sharp, i.e., 1^ <C Ir- 

Other lengths appearing in the problem can be ex- 
pressed in terms of the ones introduced above. 



III. ANALYSIS OF STEADY STATE 

Analytic treatment of the steady state is simpler if one 
neglects the hole tunneling (rh ^ oo), so that the hole de- 
pletion length /h = VDhTh disappears from the problem. 
Note that a finite Tc is necessary to provide the uniform 
electron background at large r. It is also convenient to 
first consider an infinitely sharp interface (Z„ = 0). In 
addition, the treatment in Ref. [4] neglects the electron 
density on the hole side of the interface, and vice versa. 
We will consider this simplified model in III A, first with- 
out assuming anything about In/le, and then writing out 
both In <C le and Zr » le limits. 

In IIIB, corrections due to nonzero rie in the hole side 
are derived. In IIIC, a finite Th is re-inserted, and in 
III D the width of the sharp interface itself is studied. 



A. Simplified model 



With Th ^ oo, the equations for steady state are 

DeV^He - 'JHiine +G - 



1 







-Tie = 



Te 



(2a) 
(2b) 



If the hole and electron densities are strictly zero outside 

and inside the ring respectively, then the nonlinear re- 
action terms in Eqs. (2) then drop out both inside and 
outside the ring. The resulting linear equations can be 
exactly solved: 



rih(r) 
ne{r) 



Ko {r/k) 



1 



Ko iln/le) 



r) , 

0{r - In) . 



(3) 
(4) 



Here K^{x) is a modified Bessel function of the second 
kind. 

Note that Eq. (3) is very similar to Eq. (17) of Ref. [7], 
where also a steady state is analyzed. On the other hand, 
Eq. (4) involves a length scale (/o), which is not so com- 
mon in previous studies of steady state fronts. Instead, it 
resembles more the time-dependent case of Refs. [8-11], 
where the corresponding length scale is determined by 
the time t, analogously to our Tq. 



In the Zr <C le limit, Eq. (4) can be written approxi- 
mately as 



ne(r) = Gts 



1 - 



ln(r/Ze) 
In(V^e) 



0ile - r) 



0{r - Ir) , (5) 



provided we are not interested in the ne(r) behavior for 
r > Ze- In the limit le'. 



Gt, 



Q{r - In) . (6) 



With the expressions for the densities, one can now 

match the electron infiux and hole outflux currents {j = 
—DVn) to determine the ring radius Zr in the simplified 
model: 



[ Px 1] 


= D, 


_2nDh Zr. 





Gt, Ki (Za/Ze 



Ko (ZR/Ze 



(7) 



This equation can be solved numerically to give the ring 
position Zr as a function of the hole source strength Px, 
or more "universally", to express Zr/Zc as a function of 

(Zsrc/Zc) • 

The hole diffusion constant drops out, and so the 
interface position is independent of fh- In addition, the 
recombination rate 7 does not enter becaiisc of the zero- 
width approximation for the interface. This approxima- 
tion turns out to be surprisingly good as far as Zr is 
concerned; as long as there is a well-defined peak in the 
luminescence, changing 7 affects the width and height of 
the peak profile but not the position (Fig. 2a). We also 
note that the two source parameters enter only as the 
ratio Px/G and not individually. 

In the Zr <C Ic and Zr ;» Zc limits, Eq. (7) can be solved 
analytically for Zr, giving respectively 



Zr = Ze exp [-27rGZg/Px] = Zg exp 



(Zgrc/Ze)^ 



and 



Zr = Px/27rGZe = llJ2l 



(8) 



(9) 



Eqs. (3), (5) and (8) have been obtained previously 
[4]. Ref. [4] has a spurious D^^ factor in the exponent of 
expression (8) for Zr. 

The theory developed in this section, based on the ap- 
proach of Ref. [4] , is now evaluated by comparing with 
data from numeric simulation (App. A) of the diffusion- 
reaction equations. In Fig. 2a, luminescence profiles have 
been plotted for several cases to show that the ring ra- 
dius remains unchanged if the recombination rate 7 is 
changed (an assumption of the theory), if the hole diS^u- 
sion constant £>h is changed (a prediction of the Th — > 00 
theory), and if the electron injection current density G 
and hole injection rate Px are both changed while their 
ratio Px/ G = TrZ^^j. is kept fixed (another prediction of the 
theory) . In each of these cases the width of the interface 
is modified, as discussed in Sec. HID. 
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FIG. 2; Assessment of the simple model in Sec. Ill A and in 
Ref. [4]. (a) Luminescence profiles (7nhnc) shown for situ- 
ations in which Zsrc, Ic are fixed but other parameters vary. 
Compared to the solid curve, dashed curve has 7 decreased 
by factor 25, while dotted curve has Dh decreased by a factor 
5, all other parameters remaining fixed. Dash-dot curve has 
G and Px decreased by factor 5, while keeping the ratio P^/G 
fixed, (b) Solid line is hole density nh(r < /r) in hole-rich 
side of interface, from simulation. Dash-dot line shows best 
fit of the form ~ ln(ZR/r). (c) Electron density distribution, 
nc(r- > Ib.), from simulation, normalized to and subtracted 
from nc{oo) — Gtb. Dash-dot line is fit to Ko{r/lc), Eq. (4). 
(d) Radius of ring, in units of the electron depletion length 
Ic, plotted as a function of hole source strength Px, in units 
of nGll 



In Fig. 2b, the steady state hole density profile ob- 
tained from simulation of Eqs. (1) is compared with the 
logarithmic prediction, Eq. (3). The agreement is rea- 
sonable but imperfect; an improvement will be found in 
the next section. 

In Fig. 2c, the expression (4) for the electron densities 
outside the ring, nc{r > ^l), is tested against simulation 
data. There is some deviation at large distances, which 
remains unexplained. The density profiles of Figs. 2b and 
2c are taken from a steady state solution with Zr « 0.18. 

Finally, in Fig. 2d, the radius of the circular interface, 
obtained from direct simulation of Eq. (1), is plotted 
against hole injection intensities P^. This is compared 
with Eq. (7), plotted as a solid line. The Ij^ <C 4 hmit 
is shown by a dashed line. For larger radii (Zr > and 
^src > ^o), the prediction for the radius is seriously at odds 
with the numerical results. The simulations suggest that 
the dependence on P^/G follows a lower exponent than 
the linear dependence obtained in this section. This dis- 
crepancy is corrected in the next section. 



B. Corrections from "dark" interior 

We now turn our attention to the hole-rich region 
within the ring, at radial distances r < l^, far enough 
from the reaction front so that rih 3> rig. We will now 
encounter effects of the extended source term for the 
electrons, i.e., of the position- independent G. Taking 
account of these effects turns out to be the key to over- 
coming the failure in Sec. Ill A to predict the ring radius 
for /r > Ic. 

Some of the figures pubhshed by Ref. [1, 4] suggest a 
nonzero luminescence intensity in the nominally dark re- 
gion between inner spot and ring. The small but nonzero 
intensity in the ring interior seems to be roughly constant 
between the inner spot and the ring, but a more quan- 
titative statement is hard to extract from the published 
figures. To the best of the present author's knowledge, 
this feature has not been explained previously. 

In numerical results (e.g., in Fig. 1 and also in Fig. Ic 
of Ref. [4]), one feature of the luminescence (77111 ric) curve 
is that it is nonzero and very nearly constant in most of 
the supposedly dark interior of the ring. The constant 
value is found to be equal G, the electron influx density. 
In other words, our reaction zone has an "extended" part 
in the hole side of the interface. 

To explain the constant luminescence for r < Zr, as 
seen in the numeric simulations and possibly in the ex- 
periments, we relax the assumption that ne(r) vanishes 
completely inside the ring (r < Zr). In the steady-state 
equation for the electron density, Eq. (2b) , the tunneling 
term can be neglected compared to G because no(r < 
^r) ^ nc{oo) — Gtc- The diffusion term is also small be- 
cause, away from the interface, Uc is small and smoothly 
varying. (This is justified more rigorously, a posteriori^ in 
App. B.) We are left with "fn\i{r)nc[r) w G, as required. 

A finite ^n\jic also affects the steady-state hole density 
distribution. Feeding 7nhno = G into Eq. (2a), we get a 
correction to the expression (3) for the hole density: 



nh{r < Zr) 



Px 



27rDh 



In 



G 



Px 



27rDh 



2/2 



(10) 



with ^ /r. Assuming the luminescence peak to be 
sharp enough, using the condition ni^{r — Ib.) — yields 

/J,, = /Rexp[-/|/2ec]. 

The lower inset to Fig. 3 shows that Eq. (10), with 
the (G/4I?h)''^ term included, provides perfect agree- 
ment with the numerical simulations. This can be com- 
pared to the previous attempt (Fig. 2b) . The same inset 
also shows the electron density in the hole region (much 
magnified), perfectly obeying 



nc{r < Ir) 



G/7 



nh{r<ln) ll.Hl'Jr) + ^r^ 



(11) 



Note that the decay of rio as one moves away from the 
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FIG. 3: Radius, in units of 1^, plotted against Px, in units of 
■kGII. The dots are from numerical simulations, the solid line 
is the improved theory of Sec. Ill B, and the dashed line is the 
theory of Sec. Ill A. The upper inset displays same curves 
in linear scale to give a sense of how the curves cross over 
from inverse-exponential to power-law behavior as /src crosses 
Ze. In the lower inset: Negative-slope curve is nh(r < Zr), 
fitted (thick dashed line) with Eq. (10). Positive-slope curve 
is nc{r < Zr), magnified by five orders of magnitude, with 
thick dashed line showing fit by Eq. (11). 



interface is not exponential or even power-law, but much 
weaker. 

The Z^ that gives the best fit to the nh and no curves 
is also in excellent agreement with the prediction above. 

An even more dramatic improvement occurs with the 
prediction for the radius, which we determine, as before, 
by imposing Di,\n'^{lY{)\ = Dc\n'^{lB)\: 



1 

2h 



{I: 



2 

src 



I. 



Ko{lR/le) 



(12) 



Fig. 3 shows how the peak postions obtained from direct 
simulation of the diffusion-reaction Eqs. (1) are perfectly 
explained by Eq. (12). The discrepancy of our original 
attempt following Ref. [4], as shown in Fig. 2(a), has been 
solved. 

The Zr 3> and Zr <C Ic limits are respectively 



+ (P^/ttG) = - Ze + yZl 



P 

^src 



It is remarkable that the tiny no(r < Zr), orders of 
magnitude smaller than Gt^ or nh(r < Zr), actually mod- 
ifies the global structure of the steady-state configura- 
tion. 



C. Finite hole tunneling 

We now relax the approximation of infinite hole leak- 
age timcscale Th, so that the hole depletion length Zh — 
^/DhT\ is finite and can play a role. Eq. (3) for the hole 
density is now corrected to 



nh{r) 



P. 



27rDh 



Ko 



(13) 



For r <C ^h, the Kq solution reduces to a logarithm, as 
before. 

Note that, since the Kq function does not vanish for 
finite arguments, the radius Zr cannot be built into 
nYi{r < Zr) as a boundary condition. The discontinuity in 
Eq. (13) suggests that the structure of the interface plays 
a more important role here compared to the Zh — *■ oo case. 
In addition, Eq. (13) also allows us to infer the ring radius 
Zr using "physical" arguments. The discontinuity can be 
minimized by having Zr > Zh, because the Kq{x) func- 
tion crosses over to ^ e^^/y^ for a; > 1. On the other 
hand, Zr cannot be too much larger than Zh, since the 
hole flux also decreases exponentially for Zr > Zh. The 
radius is therefore expected to be slightly larger than the 
hole depletion length Zh, for a range of parameters. 

As in Sec. IIIB, one should correct for nonzero 
7nh(r)ne(r), at r < Zr: 



nh{r < Zr) 



27rDh ° Uh 



— Grh , 



and 



nc{r < Zr) = 



211 



(14) 



(15) 



and 



Since the correction to nhij") is a constant, the extended 
part of the reaction zone loses the crucial role it had for 
Zh ^ oo in the determination of the interface position Zr. 
Assuming again an infinitely sharp interface at Zr and 
equating currents. 



Ze exp 



2Z2 



src 



P 



Ze exp 



2Z| 
P 

src 



From the solution of Eq. (12), e.g., from Fig. 3, one ob- 
serves that Zsrc/^o also tends to be large for Zr ^ Ze. 
Using this additional information, the Zr 3> Ic expres- 
sion reduces to Zr « Zsrc ~ x/Px/t^G. This explains the 
straight line in the log- log plot of Fig. 3 for large Zp/Ze. 
The line has slope half of that in the case of the simple 
theory without interior correction. Fig. 2(a), where the 
behavior is Zr cx Px- 



P 



27rZh 

In the Zr ^ Ze, Zr 3> h-i limit, 



Gl 



Ko{Ir/Io) 



(16) 



(17) 

Here Wq{x) is the principal branch of the Lambert W 
function [23]. The large-Zp behavior for comparable Zh 
and Ze is thus logarithm-like rather than power-law. 




FIG. 4: Ring radius with finite hole tunnehng. (a) Variation 
with hole depletion length Ih, for Px = 25nG. (b) Variation 
with Px/ttG, for Zh ~ V^lc- In each plot, dots are from di- 
rect simulation of Eqs. (1). Solid curves are the prediction of 
Eq. (16) and dashed curves are the Zh = cxd result, Eq. (12). 



Unsurprisingly, in the ^r, <C ^o, ^ limit, one 
recovers the ^r, ^ le hmit of Sees. Ill A and IIIB , — 
l,exp[-2Pjll^]. 

Fig. 4a shows the dependence of the radius on the hole 
tunneling. At small Zh, the radius obeys Eq. (16) well. As 
predicted, here the radius tends to be somewhat larger 
than but of the order of the hole depletion length ly^. At 
large Ih, the ring radius approaches the Zh = oo result 
of Eq. (12). There is an intermediate range of Ih where 
neither equations work. For the case shown in Fig. 4a, 
this crossover region is 21^ < Ih ^ lOlc- Presumably, an 
analytic understanding of this parameter region requires 
taking into account the interface structure details. The 
author has not been able to incorporate effects of inter- 
face structure into the prediction for the radius. 

Fig. 4b shows for — V2lc the radius as a function of 
hole source intensity. For this Ih/lc, Eq. (16) still works 
extremely well. 



D. Width and structure of interface 

The width of steady-state reaction fronts in 
diffusion-limited reaction processes is known from heuris- 
tic arguments [7, 12] to scale as {D\yDc/^JY^^ , where J 
is the flux of particles entering the interface region. 

In Fig. 5 the steady state interface width, obtained 
by simulation of Eqs. (1), is displayed as a function of 
various parameters. While the variations with the hole 
diffusion and the annihilation rate 7 do follow ±i 
exponents quite closely, the dependence on the electron 
diffusion Dq is much weaker. To understand this, one has 
to consider the particle flux J. For both cases of infinite 
and finite Th, the flux of electrons into the interface region 
is J = Glc[Ki{lYi/lc) / Kq{Iy{JIc)\- There is complicated 
dependence on the ring position Zr, but in the Zr ^ 

limit one can use [Ki{x) I Kq{x)] 1 to simplify: 

^ ^ [ iGk ) ^ jU^cy^r^/' ■ ^^^^ 
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FIG. 5: Width (full width at half-maximum) of luminescence 
peak, from steady states obtained by simulation, plotted (log- 
log) against various parameters. The error bars reflect the 
fact that one could choose half-maximum using either G or 
zero as the base value of the luminescence ^ngrih, since ynenh 
falls to G on one side of the peak and to zero on the other 
side. In (c), G is varied while keeping the ratio Px/G = tvIsic 
constant. In (e), Dh and are varied together. The straight 
lines are power-law fits with exponents (a) 1/3, for Dh, (b) 
-1/3, for 7, (c) -1/3, for G, (d) 1/6, for D^, and (d) 1/2, for 
D = Dh = Do. 



The variation of the numerically determined width with 
G (Fig. 5c) is also in accord with this prediction. In 
Figs. 5a-c the ring position Zr is unchanged. 

The variation with Dc shown in Fig. 5d is more com- 
plex; in this case Zr also changes with Dc- While the ex- 
ponent 1/6 works reasonably for an intermediate range 
of Dc, there is significant deviation at larger because 
the ring radius Zr gets smaller, leading to a breakdown 
of the Ki/Kq ~ 1 approximation. At small Dc, the in- 
terface width is difficult to define because the interface 
becomes highly asymmetric for Dg <C D^, as indicated by 
the large error bars in Fig. 5d. In Fig. 5e both diffusion 
constants are varied together. The interface width is now 
better defined over a wide range and the exponent 1/2 
(from dI/^dI^^ = 1)1/3+1/6 ^ £)i/2) ^oj.ks ^gj-y ^gU. 

The detailed structures of 7ih,c at the interface are dif- 
ficult to put in closed form. Within this region nc{r) 
crosses over from its r < Zr behavior, Eq. (11) or (15), to 
its r > Zr solution, Eq. (2b). In the same region, the hole 
density crosses over from its interior solution, Eq. (10) or 
(14), to its r > Zr solution which we have not considered 
yet. For r > Zr -f Zo, where Uc has reached ne{r = 00) = 
Gtc, the hole density decays fast, as Ko{r/lout), with the 
small decay length Zout = y/D^/GTcj = lw\/hjl^- 

We will not attempt to extract details of the crossover, 
which can in principle be obtained with a Iw/Ir expan- 
sion, similar to boundary-layer theory [24] developed in 
the context of fluid flows near boundaries. 
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IV. THEORETICAL CONTEXT AND 
LIMITATIONS 

Although motivated by particular solid-state exper- 
iments, it is instructive to consider this analysis in 
the context of theoretical investigations of steady-state 
diffusion-limited reaction fronts and closely-related situ- 
ations. A thorough study of a simple steady state front, 
with diffusion and annihilation terms and equal and op- 
posite currents, appears in Ref. [7]. Our calculations are 
in the same spirit, but we have specific source and de- 
cay terms in addition, which play crucial roles. A related 
(and more often studied) phenomenon is that of time- 
dependent fronts, where two species are initially well- 
segregated [9 11]. Many of the same considerations ap- 
ply, with powers of inverse time playing a similar 
role as the particle flux J does in the steady-state case. 
Geometries similar to ours have been considered in Refs. 
[14, 15, 25], where one species of the reaction-armiliilation 
pair forms an island in a sea of the other. 

We have limited ourselves to the mean-field diffusion- 
reaction equations (1). In principle, mean-field treat- 
ments are valid only above the critical dimension, which 
happens to be two. At and below the critical dimen- 
sion, fluctuations become important [8, 10-13]. In the 
2D system of the present Article, effects of fluctuations 
may show up in several ways. First, the form of the an- 
nihilation term we have used, "ynhrie, can be expected 
to have logarithmic corrections in 2D [8]. Logarithmic 
corrections are also expected for the power-law scaling 
of the width [8, 11, 12]. We justify the mean-field ap- 
proximation by noting that almost all the quantities we 
have considered involve large numbers of particles so that 
fluctuations are unimportant. 

We have assumed that the charged carriers annihilate 
directly, neglecting the diffusion, dissociation and quan- 
tum dynamics of bound exictons [19, 20]. Exciton diffu- 
sion might cause the observed luminescence width to be 
larger than what corresponds to jrihrie in our model, but 
the overall trends of Sec. HID are not expected to be 
affected severely. We have also ignored possible effects 
of quantum degeneracy of the charged carriers, which 
could change the form of the diffusion terms, so that 
i^h.e are themselves density-dependent [5] . All effects of 
coulomb interactions, including screening effects from the 
conducting leads [5, 6, 20], have also been left out of our 
formulation. 



V. DISCUSSION 

In their brief analytic treatment of the steady state, 
Butov et. al. [4] have assumed Zh = oo (lack of hole tun- 
neling decay) and /r <C /c- Our results of Sec. IIIC allow 
an assessment of the Ih = oo approximation (Sees. Ill A, 
IIIB and Ref. [4]). Fig. 4a shows that it is reasonable for 
III ^ lO^c, but breaks down for smaller Zh. Since Th > Te 
is typical in the experimental realizations, the Ih = oo re- 



sults may well be experimentally relevant in some cases. 

On the other hand, the Zr <C le approximation is more 
questionable. First, with the Zr ^ exp[— A/Px] behavior, 
a relatively small change in can induce an orders-of- 
magnitude change in Zr. This implies that fluctuations 
in the effective Px would cause the ring position to fluc- 
tuate wildly, so that the stable luminescence ring pattern 
would be unlikely to have been observed. (Such fluctu- 
ations have also been observed in the numerical simu- 
lations for Zr < Ze.) Second, experimental data on the 
ring radius as a function of intensity [5, 6] show power- 
law behavior rather than any strong exp[— A/Px]-like be- 
havior. While the relationship between Px and the in- 
tensity is not known, it is unlikely to compensate for 
the exp[— A/Px] behavior and give power-law-like Zr-vs- 
intensity curves. It is therefore important to consider the 
Zr Ze case in detail, as we have done. 

We now comment on the experimental Zr vs. intensity 
data [5, 6]. The non-monotonic behavior in Fig. 5 of 
Ref. [6] strongly indicates that the dependence of the Px 
parameter of our model on the laser intensity is compli- 
cated. Note that Px is an effective parameter measuring 
the amount of excess holes diffusing out of the laser ir- 
radiation region. To the best of the author's knowledge, 
the process of generating excess holes has not been mod- 
eled carefully, and nothing is known conclusively about 
the Px-intensity dependence. 

In Ref. [5], Denev et. al. have suggested that the linear 
behavior of Zr vs. intensity might be due to the impor- 
tance of coulomb terms which are not included in the 
present diffusion-reaction model. However, if the effec- 
tive Px parameter is a quadratic power of the intensity, 
our Zr oc ^/P^ prediction for Zh ^ le would also show up 
as a linear ZR-intensity result. 

Our analytic results gives insight into other simula- 
tions, for example the numerical results in Fig. lb of 
Ref. [5]. The fact that this curve behaves roughly log- 
arithmically at large Zr (large Px), rather than as a 
power law with exponent 1/2, shows that the simula- 
tions were done using finite Th, with Zh not too large 
compared to Zg. Note that the Lambert W function 
of Eq. (17) is roughly logarithmic for large arguments, 
Wo{x oo) w logx — log(log.T). 

To summarize, motivated by semiconductor lumines- 
cence experiments, we have investigated a two-species 
inhomogeneous steady state arising from mean-field 
diffusion-annihilation equations with a localized source 
for one and an extended source for the other. If the 
holes are not allowed to have single-particle (tunneling) 
decay, our analysis predicts the density profiles and the 
radius of the ring-shaped interface with spectacular suc- 
cess. When both species are allowed to tunnel out, the 
quality of the analytic predictions is more modest. We 
have detailed the crossover between finite hole tunneling 
and zero hole tunneling behaviors of the interface posi- 
tion. The thorough study of the steady state within the 
diffusion-reaction model should serve as a baseline for 
evaluating the need to invoke additional physical effects 
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for explaining experimental observations. 
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APPENDIX A: NUMERICAL SIMULATIONS 

The numeric steady states have been obtained by fol- 
lowing in time the evolution of Eqs. (1). The onc- 
dimcnsional spatial grid was not linear but chosen to 
be concentrated at smaller radial distances. The time 
evolution due to the diffusion terms was performed by a 
symmetric combination of forward and backward Euler 
evolution. This is sometimes called the "improved Euler 
method" and has error 0{6t^) per time-step. The time 
steps St themselves were determined adaptively, and kept 
small enough such that the diffusion terms would not de- 
crease densities below zero. 



The terms other than diffusion were treated "exactly" 
within each times step, i.e., to order 0{6t°°). For holes, 
the change nh{t + 5t) — nh{t) is given by 



Srih = nh(i) - -PhC 



-Ao(5t 



where Ac = "fHcit) + 1/tc acts as a decay factor. The 
electron evolution in each time step is similar with the 
source G instead of Phe~^ /'^ . 



APPENDIX B: SMALL ELECTRON DIFFUSION 
IN INTERIOR, JUSTIFIED 

To justify the neglect in Sec. IIIB of the diffusion term 
DeV^rio compared to G in the r < Zr region, one can use 
ne(r) « G/jnh{r) « 2Dh [lg,.cjln{lB,/r)] to estimate 
the diffusion term. The result can be expressed as 



l-PeV'ne 

G 



^src^R 



(^src/ 



Using the observation {Isrc/le)^ > {Ir/Is), from Fig. 3 
or Fig. 4, we see that a sufficient condition for DeW^rie/G 
to be negligible is l„ ^ Zr , which is true as long as there 
is a well-defined interface. 
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